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Abstract 

Imprinting effects can lead to parent-of-origin patterns in many complex human diseases. For hypertension, 
previous studies revealed the possible involvement of imprinted genes. Genetic Analysis Workshop 18 real data, 
with hypertensive phenotype and genotype of more than 1000 individuals from 20 pedigrees, provided us an 
opportunity to further substantiate such findings. To test for imprinting effects, we developed a pedigree-parental- 
asymmetry test taking both affected and unaffected offspring into consideration (PPATu). We carried out a 
simulation study based on the Genetic Analysis Workshop 18 pedigrees to show that PPATu has well-controlled 
type I error and is indeed more powerful than the pedigree-parental-asymmetry test (PPA"0, an existing method 
that does not utilize information from unaffected offspring. We then applied PPATu to Genetic Analysis Workshop 
18 genome-wide association study data from 20 pedigrees. We identified a number of single-nucleotide 
polymorphisms showing significant imprinting effects that are within genomic regions that have been previously 
implicated to be associated with hypertension. 



Background 

Genomic imprinting refers to the phenomenon of unequal 
expression of a heterozygous genotype depending on 
which parent (father or mother) the imprinted variant is 
inherited from. It is estimated that approximately 1% of all 
mammalian genes are imprinted (http://igc.otago.ac.nz). 
Among these imprinted genes, Beckwith- Wiedemann syn- 
drome, Silver-Russell syndrome, Angelman syndrome, and 
Prader-Willi syndrome are the best known. 

Numerous methods have been proposed to detect 
imprinting effects. For a diallelic genetic marker locus, 
the parental- asymmetry test (PAT) that considers imbal- 
ance of parental origins of the variant allele is simple and 
powerful. A series of generalizations of PAT, such as the 
pedigree-parental-asymmetry test (PPAT) for general 
pedigree data, widen its practical range [1]. However, 
these tests use only information on affected offspring and 
their parents. Recently, PATu [2] was proposed to take 
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unaffected offspring in a nuclear family into considera- 
tion, making fuller use of data to improve power. In this 
study, we propose a novel parent-of-origin effects test, 
PPATu, that uses both affected and unaffected offspring 
in general pedigrees, and apply the method to the 
Genetic Analysis Workshop 18 (GAW18) data, consisting 
of 20 large pedigrees, to study the hypertensive pheno- 
type. Previous studies revealed the possible involvement 
of imprinted genes in hypertension [3,4]. The GAW18 
data thus provide us the opportunity to further substanti- 
ate such findings. 

Methods 

Suppose that the marker of interest has 2 alleles, Mi 
and M 2 , and the disease allele is more likely to be asso- 
ciated with marker allele Mj. Let 0, 1, and 2 represent 
the marker genotypes M 2 M 2 , MiM 2 , and MiMi, respec- 
tively. For a child-parents trio, let F, M, and C denote 
the marker genotypes of the father, mother, and child, 
respectively. Throughout this article, mating symmetry 
is assumed; that is, P(F=f, M=m)=P(F=m, M=f) for all f, 
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m = 0,1,2. We also assume that there is no maternal 
effect; that is, the maternal genotype does not confer 
additional risk on the child's phenotype. 

Suppose we have N independent pedigrees, and for 
the i th pedi gree, we have nui unaffected and nai affected 

N run m 

offspring. Define s = E E fe/> M ,,c irl - Vm„c„=i) - E fe>M«,c ft=1 - h ik < Mik , c ^)] 

i=l j=l k=l 

where / is the usual indicator function. We can prove 
that under the null hypothesis of no imprinting effect, E 
(S) = 0. The unbiased estimator of the variance of S is 



i=i ;=i 



The standardized test statistic 



PPATu ■■ 



- h ih <M ik ,C ik ^)] 



follows 



the N(0, 1) distribution asymptotically. When there is 
maternal imprinting effect, PPATu will be positive; when 
there is paternal imprinting effect, it will be negative. 
Note that the contributions from trios in a pedigree are 
not independent, and their correlations are accounted for 
in the variance. In our simulation study and application 
below, we compare the performance of PPATu with 
PPAT, whose statistic is defined without the negative 
terms in the S statistic; that is, without utilizing informa- 
tion on trios with unaffected offspring. More specifically, 

N nai 

J2 J2 Oy>iWy,Q r l - ^Fy<My,Q r l) 
i=l)=l 



PPAT = 



N nai 2 

H E Oy>A4y,Q r l - ^Fy<iWy,Q r l)] 
i=l j=l 



Results 

Simulation study 

To evaluate the power of the proposed statistic and to 
compare with PPAT, we carried out a simulation study 
under 9 different settings, combinations of 3 sets of hap- 
lotype frequencies (HI, H2, H3) and 3 imprinting mod- 
els (II, 12, 13) (Table 1). Our data were simulated based 
on the general pedigrees from GAW18; their sizes are 
described below. To gauge the type I error rate, we also 
considered 9 additional settings, combinations of the 
same 3 sets of haplotype frequencies and 3 no-imprint- 
ing models (Nl, N2, N3), also given in Table 1. We 



simulated 1000 replicates under each of the settings (a 
total of 18 combinations). The results are plotted in 
Figure 1, which shows that the empirical type I error 
rates, at the 0.01 nominal significance level, are all well 
controlled for both PPAT and PPATu. On the other 
hand, PPATu is clearly more powerful in all settings, 
especially when there is a substantial imprinting effect. 
Conclusions are the same for significance levels 0.05 
and 0.005 as well; consequently,the results are not 
shown for brevity. 

GAW18 data analysis 

We consider GAW18 real genome-wide association stu- 
dies pedigree data that contain a total of 472,049 single- 
nucleotide polymorphism (SNP) genotypes on odd chro- 
mosomes and phenotype data, including systolic and dia- 
stolic blood pressure (SBP and DBP). In our study, we 
use a hypertensive binary phenotype; an individual is 
classified as affected if SBP > 140 mm Hg, or DBP > 90 
mm Hg, or on antihypertensive medication at the first 
examination. There are 20 pedigrees; the sizes range 
from 27 to 107 individuals. In total, there are only 157 
affected offspring, while there are 709 unaffected ones. 
Hence, based on the experience gained in our simulation, 
we expect a substantial gain in information for PPATu 
that makes use of information from both affected and 
unaffected individuals when compared to PPAT. To 
reduce the effect of multiple testing, we first used pedi- 
gree disequilibrium test (PDT [5]) to identify SNPs that 
are associated with hypertension at the 0.05 level, and 
then performed imprinting effect tests, focusing only on 
those SNPs. In our analysis, all trios with complete data 
within each pedigree were included in computing the test 
statistic. Furthermore, although many tests were per- 
formed, we did not attempt to correct for multiple testing 
given the small sample size (a total of only 20 pedigrees). 

Table 2 shows the cross-classification of SNPs by dif- 
ferent tests. We also provide Figure 2, which shows the 
results for a combination of significant levels. As one 
can see from Figure 2 and Table 2, PPATu detected 
more SNPs with smaller p values than PPAT. For SNPs 
detected by both methods, the results for PPATu are 
more significant. Table 3 shows the 2 SNPs (rs 12947636 
and rsl674137) that are significant at the 0.01 level for 



Table 1 Combinations of 9 imprinting settings and 9 no-imprinting settings 

Haplotype frequency 3 Imprinting model b No-imprinting model b 





DM1 


dMI 


DM2 


dM2 




<Pd/d 


<Pd/D 


(PD/d 


(PD/D 




<Pd/d 


<Pd/D 


(PD/d 


H1 


0.2 


0.0 


0.1 


0.7 


11 


0.26 


0.28 


0.37 


0.39 


N1 


0.26 


0.33 


0.39 


H2 


0.3 


0.1 


0.0 


0.6 


12 


0.24 


0.26 


0.42 


0.44 


N2 


0.24 


0.34 


0.44 


H3 


0.3 


0.0 


0.0 


0.7 


13 


0.18 


0.23 


0.53 


0.58 


N3 


0.18 


0.38 


0.58 



a Haplotype frequencies are between disease susceptibility locus(with disease allele D and normal allele d) and marker locus (alleles M1 and M2). 
imprinting models are defined in terms of penetrance probabilities: ((p d/d , (p d/D , (p D/d , <p D /d)- For the no-imprinting models,<p D/c / = cp d/D . 
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Figure 1 Type I error and power for PPATu and PPAT. The solid 
horizontal line marks the nominal significance level of 0.01. For power, 
the 9 settings are the 9 combinations of haplotype frequencies and 
imprinting models (as given in Table 1): 1 = (H1, 11), 2 = (H2, 11), 3 = 
(H3, 11), 4 = (HI, 12), 5 = (H2, 12), 6 = (H3, 12), 7 = (HI, 13), 8 = (H2, 13), 
and 9 = (H3, 13). For type I error, the 9 settings are the 9 combinations 
of haplotype frequencies and the no-imprinting models (as given in 
Table 1): 1 = (HI, N1), 2 = (H2, N1), 3 = (H3, N1), 4 = (HI, N2), 5 = (H2, 
N2), 6 = (H3, N2), 7 = (HI, N3), 8 = (H2, N3), and 9 = (H3, N3). 



Table 2 Cross-classification of results based on p values 
fromassociation(PDT) andimprinting (PPAT and PPATu) 
tests 







PPAT 






PPATu 




PDT 


<0.05 


<0.01 


>0.05 


<0.05 


<0.01 


>0.05 


<0.05 


245 


1 


12,752 


556 


59 


12,441 


<0.01 


17 


0 


798 


30 


2 


785 


>0.05 


7008 


35 


14,681 


15,007 


1328 


6682 



both PDT and PPATu. It appears that these 2 SNPs are 
novel, as they were not previously identified as asso- 
ciated with hypertension, thus further study is warranted 
to substantiate the finding. 

To gain a more global view of the extent of the role of 
imprinting in hypertension, we also carried out genome- 
wide testing for imprinting effects without restricting to 
SNPs with small p values from PDT. The analysis has 
resulted in larger sets with p values smaller than 0.05, 
which we call "significant'Tor easy reference. A search in 
the NCBI database found 49 genes implicated to be asso- 
ciated with hypertension in previous studies. Figure 3 
gives the numbers of SNPs with significant imprinting or 
association effects within these genes. We can see that 
more SNPs identified by PPATu are within these genes 
than those identified by PPAT. Specifically, 47 and 15 
SNPs within these genes are identified as having signifi- 
cant imprinting effects by PPATu and PPAT, respec- 
tively. In fact, the proportion of SNPs identified by 
PPATu that are in known genes is higher (1.5 times) 
than that identified by PPAT, although this difference is 
not substantial. Other than 1 SNP, the rest did not reach 
the threshold of significance for PDT, which may be 
explained by the power loss of association tests like PDT 
that do not account for the imprinting effect properly. 
The 1 SNP (rsl 1606492; Table 3) that yielded significant 
results from both PDT and PPATu (but not PPAT) is 
within the gene PLEKHA7 that has been implicated to be 
associated with hypertension in previous studies [6,7] . 

Conclusions 

In this article, we propose an imprinting test that uti- 
lizes both affected and unaffected individuals from 
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Figure 2 Significant imprinting results (p value < 0.05) for SNPs having small p values for association test (815 SNPs with p value < 
0.01) across all odd chromosomes. {Top) Based on PPAT that uses only affected off spring and their parent data (17 SNPs identified). {Bottom) 
Based on PPATu that uses both affected and unaffected off spring (30 SNPs identified). 
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Table 3 SIMPs having p values <0.05 for both association and imprinting tests 



SNP 


Chr 


Gene 


Position 


P-PDT 


P_PPAT 


P_PPATu 


Imprinting 


rs 12947636 


17 


SLC39A11 


70992565 


0.0090 


0.6171 


0.0080 


maternal 


rs1674137 


19 


TSKS 


50258030 


0.0089 


0.6547 


0.0038 


maternal 


rs 11606492 


11 


PLEKHA7 


17000241 


0.0488 


0.1797 


0.0176 


paternal 



The first 2 SNPs are significant at the 0.01 level for PDT and PPATu. The third SNP is significant at the 0.05 level for PDT and PPATu, and is within the PLEKHA7 
gene previously identified as associated with hypertension. 



Ppdt<0.05 



Pppat<0.05 




Pppatu<0.05 

Figure 3 Number of SNPs with significant imprinting or 
association effects (p value < 0.05) that are with in genes 
previously identified as associated with hypertension. The pink, 
yellow, and blue circles denote results from PDT, PPAT, and PPATu, 
respectively. 



general pedigrees, the type of data provided by GAW18. 
We expect PPATu to be more powerful than the existing 
test PPAT [1] because the former makes full use of infor- 
mation by taking unaffected offspring into consideration. 
Indeed, our simulation study shows that PPATu has 
higher power than PPAT without an elevated type I error 
rate based on the GAW18 pedigrees. Our results from 
analysis of the GAW18 data using PPATu leads to the 
identification of a number of SNPs that are within geno- 
mic regions previously implicated for the hypertensive 
phenotype. Nevertheless, further investigation is warranted 
especially to evaluate the performance of the methods 
under different study designs and ascertainment criteria. 
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